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Abstract 

Metals approaching the Mott insulator generate a new hierarchy in the electronic struc- 
ture accompanied by a momentum dependent electron differentiation, beyond the Mott- 
Hubbard, Brinkman-Rice and Slater pictures of the Mott transition. To consider such 
nonlinear phenomenon, we develop an analytic nonperturbative theory based on opera- 
tor projections combined with a self-consistent treatment of the low-energy excitations. 
This reproduces the formation of the Hubbard bands, Mott gap, spin fluctuations, mass 
divergence, diverging charge compressibility, and strongly renormalized flat and damped 
dispersion similar to angle-resolved photoemission data in high-T c cuprates. Main struc- 
tures in electronic spectra show a remarkable similarity to numerical results. 
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1 Introduction 

The issues of metal-insulator transitions (MIT) jlj and the high-T c cuprate superconductiv- 
ity Q have opened extensive studies on related subjects from both theoretical and experimen- 
tal aspects J|, H H D . In spite of a long history of studies || [|, (?], ^ , their appropriate and 
satisfactory theoretical descriptions remain controversial. This paper aims at giving a unified 
view of the filling-control MIT in the 2D Hubbard model by means of the operator projection 
method (OPM) || and discussing electron differentiation in the momentum space obtained by 
the OPM which also occurs in underdoped high-T c cuprates. The OPM is a recently developed 
non-perturbative analytic method based on the projeciton technique in an operator space [To[ || 
and systematically improves the Hartree-Fock theory and the Hubbard Q and two-pole approx- 
imations jyj |l2| as well as the perturbation and spin fluctuation theories |l4|, [l5|, [H| |l^] . 

Pictures of early theories on the MIT are classified into three types: (i) The Hubbard 
approximations described an MIT as the splitting of two Hubbard bands due to the strong local 
Coulomb repulsion U compared with the bare bandwidth W However, they and improved 



theories 1 1 1 , [12| violate the Luttinger sum rule in paramagnetic metals, (ii) Brinkman and 
Rice predicted that at half filling, an MIT occurs at U c (> 0) due to the disappearance of 
quasiparticles accompanied by the mass divergence However, it ignores spin correlations 
and incoherent Hubbard bands, (iii) The Slater gives a view that a folding of the Brillouin 
zone accompanied by the antiferromagnetic (AF) phase transition yields an energy gap in the 
single-particle excitations. Self-consistent renormalization theory modified the Hartree-Fock 
theory and the RPA 0. However, they do not give Mott-insulating features above the Neel 
temperature. 
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Further intensive studies on the MIT have been promoted to explain peculiar features ob- 
served in high-T c cuprates §, |, [if, ||. Their proper understanding and description remains 
still challenging due to difficulties in correctly describing strong correlations near the MIT. It 
is crucial to build a unified theory of the three aspects of the MIT, in order to explain global 
phase diagrams of compounds that undergo an MIT [jjj [|, to clarify the mechanism of the 
high-T c superconductivity and to describe the MIT in the Hubbard model. 

Gaussian fluctuations around the slave-boson mean-field theory |2(| as well as the dynamical 
mean-field solution (DMFT) Q incorporates the Hubbard bands but not spin correlations to 
the Brinkman-Rice picture. The DMFT predicts a lowest-energy resonance inside the Hubbard 
gap at half filling below U c (> 0). However, it ignores the momentum dependence of the self- 
energy part in finite dimensions. Some independent momentum patches can be inctroduced into 
this local self-energy part [|T). Then, for the 2D half-filled Hubbard model at weak couplings, 
a Mott insulator with sharp peaks at the gap edges is reproduced, though solving a cluster 
problem yields severe limitations here. 

No theories self-consistently modifying the RPA |l^, [l4|, [l5|, [ll| [lTj predict Mott-insulating 
features jl6), without abandoning the self-consistency between the Green's function and the 
self-energy part fl5|j . 

Quantum Monte-Carlo (QMC) simulations have revealed the thermodynamic and dynam- 
ical properties of the 2D Hubbard model: The charge compressibility diverges towards half 
filling Near half filling, nearly separated Hubbard bands develop, the Luttinger sum rule 
is violated at the temperature T > J = 4t 2 /U and a sharp peak of the local density of states 
develops near the Fermi level with decreasing T, reflecting a suppressed Fermi temperature 
Tf- [§3|. At T < J, extra two bands appear in addition to the Hubbard bands ^5). A 
dispersion around the (ir, 0) and (0,7r) momenta is strongly weak (^6). Strong momentum de- 
pendence of quasiparticle excitations is a crucial feature beyond the early three pictures of the 
MIT. 

The OPM |)| offers a theoretical framework unifying the early three pictures and reproducing 
key features beyond them [^|[ |23|, |25|, ^6) , the momentum-dependent electron differentiation 
in the low-energy hierarchy. They are discussed in § [| and ^|. 

The OPM is formulated in § |2|. Based on equations of motion, it gives a unique continued- 
fraction or Dyson-equation expression of a correlation function. It may systematically ex- 
tract the lower-energy properties by projecting out higher-energy degrees of freedom in a self- 
consistent fashion. For the self-energy part, we can build its Dyson equation beyond the pertur- 
bation and one-loop theories. The higher-order dynamics was considered within a two-site level 
in a local picture ]2^ | , though growing spin correlations were ignored even near the continuous 
AF phase transition. The OPM solves this difficulty Q . 

In § |^, we discuss the similarity between our results and observed normal-state properties 
in high-r c cuprates. 

§ q is devoted to the summary of this paper. 



2 Operator projection method 

We follow the previous papers ||. We consider the Hubbard Hamiltonian with an electron 
transfer t x x t(= t x , x ) from an atomic site x' to x; h = - ^ , t xx rc Xg c x , B + u^2 n xl n XL , with 
n Xs = c Xs c Xs . We perform the projection procedure for the electron creation and annihilation 
operators at a site x with a spin index s, c Xg and c Xs . It is useful to define the fcrmionic 
4jV-component vector operator \& and its Hermitian conjugate *$>\ composed of 

fx = t (ca5t,ca! 1 ,-c+. i ,ct. T ), (1) 
*i = (c xv cl r ~c xl ,c x ^, (2) 

respectively, in the real-space representation with the number of the atoms N. We define the 
commutator of an operator A with It = h - ^ n Xs as ujA = [A, H]- and the thermal average 
of A at temperature T as (A) = Ti[e^ R / T A]/ Z with the partition function Z = Tr[e~ n / T ] 
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and the chemical potential \i. The Heisenberg operator is introduced as A(t) = e Ae~~ . 
Then, for arbitrary 4iV-component vector operators with an odd number of particles and ip, 
their response function and susceptibility are defined in the AN x AN matrix representation as 
K$ >ip t(t) and X0, v t(^) = J^°dt e^K^ ^t (t), respectively. Here, the (a, a') ® (x, x 1 ) component 
of K^^t (i) is defined by 

^t(^-i([«(t),(^] + >- (3) 
The first-order projection is given by 

cjq, = e (11) * + <5cj*, (4) 
£(") = Jf^^t (0)^-^(0), (5) 
<Scj* = (1-Pxp*. (6) 

Pi is defined as the first-order projection operator which operates to an arbitrary operator 
A as Pi A = K A ! ^ t (0)K'~ 1 ^ t (0)* and satisfies P% = P 1 . <5a>* satisfies K 9 ,(s£,9)\ (0) = 
-^5(i*.*t (0) = 0. This leads to the Dyson equation for the single-particle Green's function 
G(lo) = x->& *t(w) in the 47V x AN matrix representation, 

G{u) = [(wJ-e^J-SiM] -1 , (7) 

SiH = -ixL, (M t)t(")^! tt (o). (8) 

J is the 47V x 47V identity matrix. X&» >(M ,*)tM = [x7a*,(^*)tM - i*^*t (O)G^(a;)]- 1 

is the irreducible part of XgQwJSu^i^) with respect to G^(u>). We define the local charge, 
spin and pair operators at x as n x = n xl + n xl , S x = 4 as (T aa ,c a . s ,/2, A x = -^={c xl c xl H-c^c^) 
and A x = -^(c xr c xl - c^c^). We introduce the identity and Pauli matrices uo {po) and <r (p), 
respectively, which operate to the spin (charge) space of vector operators. Then, for U > 0, 
one obtains 

= ^L'^ero - t/(Si><5 aJ>aJ /po^3, (9) 

4°,x' = -w-tc-^w/ 2 )^. ( 10 ) 

<5w*a! = f ( y Sn x p 3 <j„/2 - <5Sjp CT 3 ) *as, (11) 

with 5n x = n x - (nx) and 5SJ, = SJ - (s x ). We have taken (A^.) = o, and also (s x ) = o without loss 
of generality. 

The present theory without Si(w) yields the Hartree-Fock theory. Self-consistent calcula- 
tions of by perturbation or one-loop theories |0|, [l4j do not reproduce the Mott insula- 

tor mm. 

Next, we construct the Dyson equation for Si (a;), performing the second-order projection; 

ljSuj^ = £ (21) * + £ (22) (5w* + (5wfc*, (12) 
E ( !1 )EK iK „ t (0)^ t (0), (13) 

£ ( 22 ) = K QSO y >(smV (0)K-^ {s ^ y . (0), (14) 
SuiSluV = (1 - P 2 )Cj5uj^. (15) 

P2 is defined as the second-order projection operator which operates to an arbitrary operator 
A as P 2 A = PiA + K A ^ SSl ^(0)Kj} 9 ((5( - w)t (0)5w* and satisfies P 2 2 = P 2 . Then, we obtain 

Si (a;) = [( w 7- £ ( 22 ))-S 2 H]- 1 £( 21 ), (16) 
s 2(^) = -i7dsasu9,(6u,s^( u ) K 7ii%(s^( )- ( 17 ) 

Here ' &a*,(j«a««a*)tH = \x^66>v,(,ss>se>9)i( u ) -^^(^tW^M! -1 is the irreducible 
part of /t5w«w*,(5ii5ii*)t(w) with respect to G^(o;) and Si(w). For the repulsive Hubbard 
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model, 



= U'^M*, (18) 
Mm = {(Sn xP3 a /2 - SS^po^s) 2 ) (19) 

e a!,x' = ~~ t x%' P3<70 ~ (^ 2 ^ 3<T o - (S^poo-s) S x a .t, (20) 

t {2 % =t x ^,(8C xP3 8C x ,p 3 )M x }, (21) 

iCa; = <5n x p 3 o- /2 - (5.Sa; ■ p (T + A^pjao/V^, (22) 

= pf' pa<?a + [(1 - <n))£ki„ + £cor2] M~ 1 , (23) 

£ kin = ^■ t fe((4/fc^ P3,T0 " < C L-' T355 ' C fe5'>' ,C " T3 )' ( 24 ) 

where p^ 0) = p - c/(i - (n x )/2). Hereafter, the summations over i = 1 and 2 and barred variables 
as x and fe should be taken. Below, a local energy shift e CO r2 due to two-site correlated hopping 
terms is neglected. We note that in the particle-hole symmetric case, £ co r2 vanishes. Then, in 
the case of (S x ) - o, we obtain 

S&S&iVa W («a.,a(l - <«))e kin M x 1 + i^) p 3 5(D* a 

- Utv&SCaspaVm- (25) 

Without S2(w) and the momentum dependence of e^ 22 \ this theory gives the Hubbard I 
approximation pj. The momentum dependence of e^ 22 ) that includes a superexchange term 



modifies the dispersions of the Hubbard bands 12 



In the symmetry-unbroken phase, we have only to discuss the spin-independent electron 
part, E 2 o(w,fe) = 4XsaYa 9t (s&s&v)f(u,k)/(ri)(2 - (n)). We take the following decoupling ap- 
proximation to E2 C (w, fe) to correctly estimate its moment and to include the short-ranged 
correlations neglected in previous equation-of-motion approaches j|, [12, p7| ; 



E 2e (w„,fe) = - JJ N V [G(w n -fi m ,fe-Q) 
(n)(2-(n)) f-^ L 



(Xc(a„, <z) + Xs(On, 9)) + 24_ qX s(« m , g)} 
-2tS_,G(fi ro -w n , g-fe)x P (O m , g)] , (26) 

where T fc;(J = t fe _ g - ^J^ e^ - t k with B kin = ^ £ fe s t fc < c t, s c fcs ) and t fc is the (1, 1) component 
of t^ 22) . f2 m = 2ttitiT is a bosonic Matsubara frequency, G c is the spin-independent electron 
Green's function and Xc Xs and Xp are the charge, spin, local-pair susceptibilities, respectively. 

We calculate the susceptibilities with the two-particle self-consistent method |l5[ to include 
short-ranged correlations and to reproduce their local moments. Following Ref. (lq j, we take 

Xc(a«, fe) = 2x P h(ft m , fe)/[l+r c x ph (O ro , fe)], (27) 

Xs (n m , fe)=2 Xp h(n m , fe)/[i-r sXp h(n m , fe)], (28) 
x P (n m , fe)=x P p(o m , fe)/[i+r P x PP (o m) fe)], (29) 

with an approximation r s = iu(n xr n xl ) /(n) 2 and the self-consistency conditions for T s , T c 
and T p given by the local moment sum rules; k x»,»(^m> fe) = <™> + 2(nan n au) - (n) 2 , 

T/wJ] mife Xsi,si( n "«'fe) = <") -^i X} n xl ), and r/iV^ m _ fe x A , A t( n m.fe) = (nan«au>- Xph and x PP are 
the particle-hole and particle-particle susceptibilities calculated from the bare Green's functions. 

Below, results are shown for the single-particle density of states p(to) = ^ X)fc ^( w i fe)> 
spectral functions A(u,k) = — -llmG (cL;, fe), dispersions and the fi-(n) curve. 



3 Mott-insulating properties at half filling 

Figure [j] shows that at low temperatures, a Mott-insulating gap grows, indicating the Mott- 
insulating ground states. This is accompanied by sharp peaks of p(to) at the low-energy edges 
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Figure 1: The local density of states p(ui) at half filling (n) = 1. The Mott-Hubbard pseudogap 
remains even at T = J ee 4t 2 /U above which AFM spin correlations are not enhanced. 




Figure 2: A(cu, k) for U/t — 4 in the case of half filling (n) = 1. 



of the Hubbard bands. This agrees with the QMC results p3| but contrasts with the Hubbard 
approximations j| and the DMFT js) . As shown in Figs. | and 0, the growth of AF spin 
fluctuations develops the AF shadows of the Hubbard bands, yielding a four-band structure, 
in agreement with QMC results |2q] . The structure takes the form of a superposition of two 
SDW-like bands and two low-energy narrow bands at the gap edges. The low-energy bands 
have a particularly weak dispersion and dominant weights around (tt, 0) and (0, ir). This is 
consistent with the QMC results ^6|. These excitation spectra around (%, 0) and (0, 7r) 

mainly contribute to these sharp peaks. With increasing T, spectral weights fill in the gap and 
the intensities of the peaks are weakened. AF spin correlations also diminishes. This smears 
out the AF shadows. Then, the main structure is characterized by the Hubbard bands nearly 
separated by the Mott-Hubbard pseudogap, as shown in Fig. ||. 

For the momentum distribution, the error of the OPM compared with the T = QMC 
results turns out to be less than 10% for U/t = 4 and is even smaller for U/t = 8 p6fl. 



Moreover, Fig, p 
QMC results [p6 



suggests that a main structure in A(uj, k) shows a remarkable similarity to the 
. This small discrepancy mainly originates from a small reduction of the Mott 
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Figure 3: Single-particle dispersions for U/t = 4 and T/t = 0.02 at (n) = 1. The squares and 
the red circles denote the present and QMC results respectively. The open and the filled 
symbols represent the momenta where the peak intensities divided by the largest intensity are 
less and larger than 0.1, respectively. The solid, dotted and dashed lines denote the free-electron 
dispersion and the upper Hubbard dispersion by the Hubbard I approximation and the two-pole 
approximation with t k = t k+ Q, respectively. 



gap from U compared with QMC data. 



4 Metallic properties in doped cases 

Results on p(u>) at (n) = 0.87 are shown in Fig. ^. The upper and lower Hubbard bands are 
nearly separated by the Mott-Hubbard pseudogap with remaining weak incoherent spectra. This 
is also shown in Fig. |5| that gives A(w, k) similar to the results of the two-pole approximation [12]. 
With decreasing T, a sharp peak of p(uj) develops near u = 0. They agree with QMC results [23]. 
As shown in Figs. ^| and ^, the single-particle dispersion in the lower Hubbard band is strongly 
weak both outside the magnetic Brillouin zone boundary and around the (tt, 0) and (0, tt) 
momenta at T/t = 0.02, in agreement with QMC results ^3|. Particularly, Fig. [| shows a 
formation of a low-energy narrow band pinned at the top of the lower Hubbard band. This 
band has an extremely weak dispersion and dominant weights around (jr, 0) and (0, tt) and 
seems to evolve into the low-energy band obtained at half filling. This contributes to the sharp 
peak of the density of states at the top of the lower Hubbard band which persists even in the 
Mott insulator. This generates a divergence of the compressibility as (n) — > 1, as shown in 
Fig This contrasts with the Brinkman-Rice picture where the compressibility remains finite 
since the sharp peaks are absent in the Mott insulator ||. 

The main peak at every momentum on the magnetic Brillouin zone boundary is located 
below the Fermi level for U/t = 8 and (n) = 0.87 even at T/t — 0.02. This gives a larger 
Fermi volume than the Luttinger volume. Results on the momentum distributions also reveal 
these features ||. Here, we note that the obtained Fermi volume decreases as T — > 0. We also 
note that the present method gives a Fermi volume consistent with the Luttinger volume for 
U/t = 1 and (n) = 0.87 at T/t = 0.02. Suppressions of the Fermi degeneracy and thus Tp due 
to strong local repulsion in the proximity of the Mott insulator may produce the inconsistency 
with the Luttinger theorem for U/t — 8 and (n) — 0.87. Actually, QMC calculations show that 
for T/t ~ 0.5, the (tt, 0) level is below the Fermi level The strongly flat dispersion at low 
energies makes it difficult to speculate the T — properties from the present results. 

Figure @ shows the mass divergence around both momenta (tt, 0) and (7r/2, tt/2), although 
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Figure 4: The T dependence of the local density of states p(uS) at n = 0.87 for U/t = 8. The 
inset shows the U dependence at the same filling n = 0.87 at T/t = 0.02. A sharp peak develops 
near uj = at low T. 



the anisotropy is large. The diverging mass itself is consistent with the Brinkman-Rice pic- 
ture ||, reflecting the generic suppression of the Fermi degeneracy near the MIT. However, 
the enhanced density of states due to the low-energy band is pinned at the top of the lower 
Hubbard band. This is the origin of new physics not contained in the Brinkman-Rice picture. 
Furthermore, damping rates 7 fcp exhibit a prominent increase as (n) — ► 1. Therefore, towards 
the MIT, quasiparticles become ill-defined at least at T = 0.02i, since 7 fep ~ v F . 

We summarize a unified view of the MIT in the 2D Hubbard model beyond the Mott- 
Hubbard, Brinkman-Rice and Slater pictures: 

Towards the MIT, the local spin moments develops. It increases with increasing U through 
the reduction of the double occupancy. This produces the splitting of the upper and lower 
Hubbard bands as in the Mott-Hubbard picture H, Q]. At T < J, the large local moments 
tend to align antiferromagnetically to gain a superexchange energy, as in the Slater picture. 
This contributes to formation of the SDW-like dispersions in the Hubbard bands through the 
momentum dependence of the self-energy part. At low temperatures, a narrow band appears at 
the top of the lower Hubbard band for hole-doped cases. However, the origin of this narrow band 
contrasts with the Brinkman-Rice theory jf| and the DMFT, because the enhanced density of 
states already exists at the gap edge in the Mott insulator. Since the dispersions in the Hubbard 
bands take the form of SDW-like dispersions, the level of the narrow band is located near the 
magnetic Brillouin zone in contrast with the Brinkman-Rice picture. The velocity of these 
low-energy excitations increases with increasing holes. But the narrow band is pinned at the 
top of the lower Hubbard band, in contrast with the results of the DMFT. These lead to the 
diverging charge compressibility. 



5 Comparison with high-T c cuprates 

In Bi2Sr2CaCu208+,5, a strongly weak dispersion both around the (tt, 0) and (0, 7r) and along 
(n, 0)-(tt, 7r) has been observed with the angle-resolved photoemission spectroscopy (ARPES) |19|, 
PS| . The reduction from the band-calculation data is significant. This flatness depends only 
weakly on both T and the hole concentration; it persists below and above T c . This suppression 
of the dispersion is attributed to strong correlation effects related to the MIT, as discussed 



previously |23|, 26 1. The OPM gives an intuitive explanation of such enhancement of the elec- 
tron mass: Local AF or singlet correlations increase \t x x ,\ between the nearest-neighbor sets 
of x and x' . This affects the dispersion in the level of the two-pole approximation so that the 
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Figure 5: A(ui, k) for U/t = 8 and (n) = 0.87. 




Figure 6: Single-particle dispersions for U/t = 8 and (n) = 0.87 at T/t = 0.02 with the same 
symbols as in Fig. |L 

(n, 7r) level in the lower Hubbard band shifts downwards and the dispersion around the (tt, 0) 
and (0, 7r) momenta becomes natter. Previous analytic theories do not reproduce such strong 
correlation effects. The similarity of the dispersions between our results and high-T c cuprates 
indicates that to capture essential features of the superconductivity it is crucial to describe the 
correlation effects with strongly momentum-dependent renormalization shared by the Hubbard 
model and the cuprates near the MIT. We speculate that the low-energy narrow band plays an 
important role in recovering the quantum-mechanical coherence. 

6 Summary 

The 2D Hubbard model at and near half filling has been studied by the OPM. This is a non- 
perturbative analytic theoretical framework improving the Hartree-Fock theory, the perturba- 
tion and one-loop theories of the self-energy part and the Hubbard and two-pole approxima- 
tions |^]. The OPM captures crucial ingredients of the MIT derived from the Mott-Hubbard, 
Brinkman-Rice and Slater pictures together with the strong momentum dependence of the 
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Figure 7: fj, — U/2 versus the hole concentration 1 — (n). Near half filling, k = d(n)/dfi appears 
to diverge as (n) — > 1. For U/t = 4 and 8, the chemical potential shows a strong (n) dependence 
for (n) > 0.95 and 0.9, respectively. This occurs because slight but finite spectral weights still 
remain inside the T = Mott gap of order of U due to thermal effects. Because of the smaller 
residual states for larger U, /i dependence of (n) increases with increasing U. 

single-particle spectra, which emerges as a crucial aspect of the strong correlation and a source 
of electron differentiation in the momentum space nemurically obtained and also observed in 
high-T c cuprates. Further studies on the low-energy hierarchy are necessary to discuss the 
Fermi-liquid properties and the mechanism of the high-T c cuprate superconductivity. 

We thank F. F. Assaad for valuable discussions and providing results of Monte-Carlo cal- 
culations. The work was supported by the Japan Society for the Promotion of Science under 
grant number JSPS-RFTF97P01103. 
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